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ABSTRACT 

Aims. We investigate the molecular and dusty environment of OH/IR stars in order to characterize the mass-loss process during the 
tip-AGB superwind phase. 

Methods. Employing the AMBER instrument at the VLT Interferometer we obtained near-infrared H- and A'-band spectro- 
interferometric observations of the three OH/IR stars IRAS 13479-5436, IRAS 14086-6907 and IRAS 17020-5254 with a spectral 
resolution of about 35. We use a two-component geometrical model, consisting of a uniform disk and a Gaussian disk, to obtain 
characteristic angular sizes of the central stellar sources and their dust envelopes, as well as the flux ratios between these components. 
Results. Angular uniform disk diameters of the three central components of the objects above have values between 3.2 mas and 
5.4 mas. For their dust envelopes, we find FWHM values between 17.1 mas and 25.2 mas. The three objects show significantly dif- 
ferent flux contributions of the shells to the total near-IR flux of 61%, 38%, and 16% for IRAS 13479-5436, IRAS 14086-6907, and 
IRAS 17020-5254, respectively. According to distance estimates from the literature, the central stellar components have radii between 
900 Ro and 1400 Rq, while their dust envelopes reach FWHM values between 9000 Rq and 13000 Rq. The visibility functions of all 
three sources exhibit wavelength variations that resemble those of earlier VLTI/AMBER observations of semi-regular and Mira vari- 
able AGB stars. These are interpreted as characteristic of atmospheric molecular layers lying above the photosphere. 
Conclusions. The derived characteristic sizes of both, the central stellar atmospheres and dust envelopes are consistent with the 
canonical properties of OH/IR stars. The spectral visibility variations resemble those of other AGB stars and indicate the presence 
of molecular layers, confirming that these are a common phenomenon among AGB stars of very different luminosities and mass-loss 
rates, alike. We also find that the dust envelopes have a clearly larger optical depth than those known for Mira stars. We interpret this as 
an expected result of the "superwind" phase, the final 10 000 to 30 000 years of AGB-evolution, when the mass-loss rate increases by a 
factor of 10-100. By their different optical depths, the three dust shells studied here may represent different stages of the "superwind" 
and different initial masses. 
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, 1. Introduction 

. OH/IR stars are medium-mass, highly evolved objects, which ei- 
■ ther have nearly reached the tip of the asymptotic giant branch 
(AGB) or may just have left the AGB toward the post-AGB 
stage. Observational characteristics are an infrared excess and 
the presence of hydroxyl (OH) masers, as well as variability due 
, to the release of cool, dust-driven winds with a high mass-loss 
rate ("superwind"), prior to forming a planetary nebula (PN). 
These critical, late phases are of fundamental interest to a com- 
plete understanding of stellar evolution. 

Mass-loss rates of stars at the tip of the AGB c an reach val- 
ues of up to ~ IQ-'^Mq yr 1 (iDe Beck et al.ll2010l) . Such a "su- 
perwind" can expel a considerable fraction of the stellar mass in 
a relatively short tim e (< 100 OOP years). Time-dependent, dust- 
driven wind models dWachter et al.ll2002t iMattsson et"al]i2010l) 
reproduce such mass-loss rates reasonably well for carbon-rich 
AGB stars. A pulsating photosphere for the initial mechanical 



* Based on observations made with the ESO VLT Interferometer at 
the La Silla Paranal Observatory under programme ID 081.D-0325. 



energy input and low effective temperatures lead to conditions 
favorable for dust formation in the denser regions of such a cir- 
cumstellar envelope (CSE). Due to their high absorption coeffi- 
cient, carbon dust grains are accelerated by the stellar radiation 
such that the star is eventually surrounded by an expanding en- 
velope. It significantly changes the spectral appearance of the 
central star, due to its large optical depth at visual and near-IR 
wavelengths and thermal re-emission in the farther IR. Except 
for the case of the dust-enshrouded Mira variable IRC -20197 
(iJeong et al.ll2003l) . such self-consistent wind models do not yet 
exist for oxyg en-rich AG B stars, and, moreover, re cent stud- 
ies bv Woitk^^ ( 120061) and lHofner & AndersenI (l2007h indicated 
fundamental problems in our understanding of the wind-driving 
mechanism for oxygen-rich AGB stars due to an insufficient ra- 
diation pressure on silicate grains. 

However, previous high spatial resolution work on OH/IR 
stars confirmed the presence of circumstellar dust shells of large 
optical depth for oxygen-rich stars near the tip of the AGB. K- 
band speckle interfer ometry of the OH/IR star OH 104.9H-2.4 
dRiechers et alJ l2004l) measured a circumstellar dust shell of 
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Table 1. Characteristics of the sources 



IRAS Name 


K 


J-K 


K-L 


D 


S OH 












(kpc) 


(Jy) 


(km s-') 


13479-5436 


4.89 


4.39 


2.15 


2.9 


0.68 


17.7 


14086-6907 


4.27 


3.28 


1.98 


2.4 


6.15 


13.4 


17020-5254 


4.15 


2.59 


1.38 


2.5 


1.08 


13.4 



large optical depth of 6.5 at 2.2 yum with an apparent angular di- 
anieter of 47 + 3 mas (FWHM), corresponding to 112 + 13 AU. 
IChesneau et al.l (l2005h resolved the mid-IR dusty envelope of the 
OH/IR star OH 26.5-1-0.6 to a FWHM size of 286 mas x 214 mas, 
or about 390 AU x 290 AU. Size and complexity of the circum- 
stellar structures of OH/IR stars narrate the mass-loss history of 
the final AGB stages, which in turn shapes the future planetary 
nebula. The recent development of very high angular resolution 
instrumentation has finally enabled us to decipher this crucial 
evidence. 

In this research note, we study the three OH/IR stars IRAS 
13479-5436, IRAS 14086-6907 and IRAS 17020-5254. All 
three objects exhibit a typical two-peak maser profile at 1612 
MHz ([teLintel Hekkert et al. 1991), which is indicative of an ex- 
panding CSE. Furthermore, a large infrared excess and a strong 
9.7/im silicate emission feature in the mid-IR spectra indicate 
the presence of an optically thick, dust-ric h CSE. Table |T] show s 
the near-infrared photometry and colors of Fouque et alT dl992h . 
the bolometric distances given by Lepine et al. ( 1995), the aver- 
age flux densities Son at 1612 MHz between the two OH maser 
peaks, and the res pective expansion velocities Vexp of the OH 
shells according to lte Lintel Hekkert et a n (ll991h . 



2. Observations and Data Reduction 

We obtained near-infrared H- and /T-band spectro- 
interferometric observations of the three OH/IR stars IRAS 
13479-5436, IRAS 14086-6907 , and IRAS 170 20-5254, us- 
ing the AMBER instrument (Petro vet all l2007h of the VLT 
Interferometer (VLTI) in low resolution mode (R ~ 35) and 
three of the 1.8 m Auxiliary Telescopes (AT). These were posi- 
tioned on stations EO, GO and HO, providing ground baselines 
of 16 m, 32 m and 48 m, along the same ground position angle 
of -109° East of North. The AMBER low resolution data also 
include coverage of the J band, but because of low J-band fluxes 
and, hence the poor quality of the respective visibility data, we 
were unable to include the J-band in our study. 

Because our target stars are dust-enshrouded objects, they 
have V magnitudes too faint to use them as Coude guide stars. 
Instead, we resorted to the off-axis guide-star option of the 
VLTI and chose the stars GSC 0866800251, GSC 0924401731, 
and GSC 087 2601051 from the HST Guide Star Catalog II 
(iLasker et al.l f2008). These guide stars are located at distances of 
42 arcsec, 39 arcsec, and 55 arcsec, respectively, from our target 
stars listed above. The external VLTI fringe tracker FINITO was 
not used. The AMBER detector integration time (DIT) was set at 
0.1 sec. After each target a calibration star of well known diam- 
eter from the catalog of Merand et al. (2005) was observed (HD 
126284 with ©UD = 1-14 + 0.02 mas, HD 114912 with 0ud = 
0.95 + 0.01 mas, and HD 154486 with ©ud = L06 + 0.02 mas). 
The observations presented here were all obtained on 4 April 
2008. Table 12] summarizes the observational details. 

Data reduction was performed with standard procedures of- 
fered by the amdlib 3.0.3 package, as provided by the AMBER 
consortium and the Jean-Marie Mariotti Center dTatuUi et alJ 
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Fig. 1. Squared visibility amplitudes as a function of wavelength 
of IRAS 13479-5436, IRAS 14086-6907, and IRAS 17020-5254 
(from top to bottom). The position of H2O and CO bands are in- 
dicated. The blue dashed lines show the best-fit of a single com- 
ponent uniform disk model. The red solid lines indicate the best- 
fit two-component model consisting of a uniform disk describing 
the central source and a Gaussian disk describing the circumstel- 
lar dust shell. The fit parameters of the two-component model 
are indicated in Tab. [3] In each panel, the three sets of lines cor- 
respond to baselines (from top to bottom) EO-GO, GO-HO, EO-HO 
with increasing projected baseline lengths as listed in Tab. |2] 



120071: IChelli et al.ll2009l) . The wavelength scale was corrected 
by a scaling-factor of 1 .03 and a linear offset of -0. 14yum. These 
parameters were obtained by a comparison of the observed flux 
distri bution of the calibrators with a synthetic flux curve as in 
Wit tkowski et aP (1201 11) . The synthetic flux curve includes the 
telluric spectrum, the stellar spectrum, and transmission curves 
of optical components of AMBER. Visibility spectra were ob- 
tained by using an average of the two transfer function mea- 
surements, which were taken closest in time to the science tar- 
get measurements. Fig.[T]shows the resulting visibility data as a 
function of wavelength and Fig.|2]shows the same visibility data 
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Table 2. Observation log of 4 April 2008 



Target 


R.A 


Dec 


Starting time 


Category 


Projected baseline 


RA. 


Seeing 












(m) 


(deg) 


(arcsec) 


IRAS 13479-5436 


13 51 12.6 


-54 51 13 


05:15:18.650 


Science 


15.89/31.76/47.65 


-110 


0.83 


HD 126284 


14 26 22.2 


-55 15 33 


05:39:33.025 


Calibrator 






1.01 


IRAS 14086-6907 


14 12 50.5 


-69 21 10 


06:05:21.233 


Science 


15.49/30.97/46.47 


-100 


0.86 


HD 114912 


13 15 25.5 


-69 40 45 


06:33:44.509 


Calibrator 






1.47 


IRAS 17020-5254 


17 06 00.7 


-52 58 47 


06:58:34.476 


Science 


15.92/31.83/47.76 


-128 


0.96 


HD 154486 


17 08 08.2 


-48 53 01 


07:24:43.020 


Calibrator 






0.95 



Table 3. Fit results of the two-component model 



IRAS 


13479-5436 


14086-6907 


17020-5254 


UD diam. [mas] 


4.2 ± 


0.6 


5.4 ± 


0.4 


3.2 ± 0.4 


CSE, FWHM [mas] 


18.1 d 


:0.9 


25.2 4 


:4.4 


17.1 ± 3.3 


%flux ofUD 


0.39 ± 


0.02 


0.62 + 


0.03 


0.84 + 0.02 


% flux of CSE 


0.61 ± 


0.02 


0.38 + 


0.03 


0.16 + 0.02 



as a function of spatial frequency, together with visibility models 
discussed below in Sect. [3] 



3. Modeling the Visibility Function 

As a first approach to describe the visibility data, we used a fit of 
a single uniform disk (UD). A single-component model would 
be expected for either a stellar source that does not show a sig- 
nificant near-IR CSE, or for an optically thick CSE that com- 
pletely obscures the central stellar source. The formally best- 
fitting curves are indicated in Figs. [T]& |2]by the dashed blue 
lines. Compared to the visibility data, these curves clearly illus- 
trate that a single-component model can not provide a good de- 
scription of our sources. In particular, a single-component model 
can not explain the comparable small differences of the visibility 
data for different baseline lengths. 

We then used a simple geometrical two-component model 
of a UD representing the central stellar source and a Gaussian 
disk representing the dusty CSE. The free parameters of this 
two-component model are the angular diameter of the UD, 
the FWHM of the Gaussian disk, and the flux ratio be- 
tween these two compon ents. Following iDvckl ( l2000l) and 
iBerger & SegransanI (l2007h . let 6\jo be the angular diameter of 
the central UD, ^oauss the FWHM of the Gaussian disk, and Vp 
the fraction of the flux contributed by the UD, i.e. the fraction 
of the attenuated stellar flux shining through the dusty envelope. 
The visibility amplitude V{s) as a function of spatial frequency 
s is then given by: 



V(s) = Vp 



+ (1 - Vp) 



exp - 



(1) 



where Ji is the first-order Bessel function of the first kind. A sat- 
isfactory fit to the visibilities was obtained using a Levenberg- 
Marquardt least-squares minimization algorithm. The best-fit 
parameters of the two-component model are listed in Table [3] 
Because we use all spectral channels in the fit procedure, the 
resulting values are near-IR values averaged over the H- and 
/T-bands. The best-fit curves are indicated by the red solid lines 
in Figs. [T] & 12] These curves indicate that the two-component 
model provides a good description of the overall shape of the 
visibility curves of all three sources. 

In comparison to the two-component model, the measured 
visibility data of all three sources show additional wavelength- 
dependent features that are not explained by our simple UD 



model with wavelength-independent diameter. The visibility val- 
ues show local maxima near wavelengths of 2.2yum in the K- 
band and near 1 .7 yum in the //-band, and decreasing visibil- 
ity values toward the edges of the H and K bands. The de- 
crease of the visibility values coincides with the positions of 
H2O and CO molecular bands, as indicated in Fig. [T] These 
visibility features resemble those that have previously been de- 
tected in AMBER data of s everal Mira variable AGB stars 
(IWitfl<:owski e t al. 2Q0 8ll2011h , of the super-giant or super- AGB 
star VX Sgr (.Chiavassa et al. | 20101). and of the semi -regular 
AGB star RS Cap (Martf-Vi dal et aLlboi ih . These visibility fea- 
tures were interpreted as being indicative of molecular layers, 
in the near-IR most importantly H2O and CO, lying above the 
continuum-forming photosphere. In other words, this means that 
our simple UD model of the central stellar source could further 
be sub-divided into a wavelength-independent photospheric con- 
tinuum diameter and a close atmospheric molecular shell with 
wavelength-dependent molecular opacity. The presence of the 
same features in the AMBER data of all three OH/IR stars of 
our sample confirms that atmospheric molecular layers are a 
common phenomenon among AGB stars of very different lumi- 
nosities and mass-loss rates, alike. We note that the best-fit UD 
diameters listed in Tab. [3] are then expected to overestimate the 
continuum photospheric radii, as a result of the additional inten- 
sity contribution by the atmospheric molecular layers at larger 
radii. For less evolved AGB stars, this overestimate niay ajnqunt 
to up to about 30%, depending on the stellar phase (lire land et'al] 
2004). 

Ideally, of course, we should match the visibility data with 
synthetic visibilities derived from physical models of OH/IR 
stars. However, consistent dynamic model atmosphere and wind 
models are not yet av ailable for oxygen- rich AGB stars, ex- 
cept for the attempt by Ijeong et aP (l2003l) . Also, dust-free dy- 
namic model atmospheres do not exist for the stellar parameters 
of OH/IR stars exhibiting luminosities above lOOOOL© and ef- 
fective temperatures below 2000K (Lepine et al. 1995), so that 
approaches of a combination of dust-free dynamic model at- 
mospheres with radiative transfer models of the dust shell can 
not be attempted (as performed for Mira variab le AGB stars by 
Wittkow ski et al ] |2007I: iKarovicova et al.ll20Tll) . This means as 
well that our fit results of the two-component model in Tab. [3] 
are only indicative of the general characteristic dimensions of 
the central stellar component and the dust shell, but do not pro- 
vide accurate radius determinations of certain layers. 



4. Results 

We here summarize the best-fit parameters of the above- 
mentioned two-component model. In order to estimate absolute 
sizes of the central stellar component, including the continuum- 
forming layers and overlying molecular layers, and the outer cir- 
cumstellar dust shells, we rely on the bolometric distances given 
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Fig. 2. Squared visibility amplitudes as a function of wavelength 
of IRAS 13479-5436, IRAS 14086-6907, and IRAS 17020-5254 
(from top to bottom). The blue dashed lines show the best-fit of 
a single component uniform disk model. The red solid lines in- 
dicate the best-fit two-component model consisting of a uniform 
disk describing the central source and a Gaussian disk describ- 
ing the circumstellar dust shell. The fit parameters of the two- 
component model are indicated in Tab. |3] 



by iLepine et al.l (1 19951) (see Table [U, albeit their uncertainties 
may amount to up to a factor of two. 

IRAS 13479-5436: The best-fit visibility model consists of a 
CSE with Gaussian FWHM of 18.1 + 0.9 mas plus a central UD 
of 4.2 + 0.6 mas in diameter. The UD contributes 39% of the ob- 
served near-IR flux of this object. The angular scales correspond 
to a near-IR FWHM of the CSE of ~ 1 1 300 Rq or -53 AU and a 
radius of the central stellar disk of ~ OOOR©. For comparison, 
th e calibration of the s tellar radius as a function of K - L color 
bv 'Lepi ne et al.l ([1995) gives a stellar radius of ~ 17OOR0. The 
rather small contribution from the stellar UD disk, as well as the 
infrared colors (largest J - K color in Table [TJ, indicate a large 
optical depth of the CSE. 



IRAS 14086-6907: The best-fit visibility model consists of 
a CSE with Gaussian FWHM of 25.2 + 4.4 mas plus a central 
UD of 5.4 ± 0.4 mas in diameter, which contributes 62% of the 
total flux. The corresponding FWHM of the CSE is -BOOORq 
or ~60AU, and the radius of the central UD is ~14OOR0. For 
comparison, the K - L calibration by Lepine et al. gives a stellar 
radius of ~1500Ro. In terms of both, stellar flux contribution 
and J - K color, this object consistently takes the place between 
the other two. 

IRAS 17020-5254: The best-fit visibiHty model consists of 
a CSE with Gaussian FWHM of 17.1 + 3.3 mas plus a central 
UD of 3.2 ± 0.4 mas in diameter, contributing 84% of the to- 
tal flux. The corresponding FWHM of the CSE is -9000 Rq or 
~43 AU, and the radius of the central UD is ~86OR0 (compared 
to ~912Rq derived from the K - L calibration). The large flux 
contribution from the stellar disk, as well as the small value of 
J - K - 2A, both indicate that we see mostly the stellar disk in- 
cluding atmospheric molecular layers, with a comparably small 
contribution from a diluted, fairly transparent dust-envelope. 

Please note that there is a consistent and, hence, convincing 
relation between (larger) IR-colors and (larger) flux contribution 
from the dust-rich CSE of these three objects. 

5. Discussion 

From computati ons of hydrodynami cal wind models of long- 
period variables dWinters et alj|2000h we derived CSE radii, as 
they would appear in the near IR, to compare to the observed val- 
ues. Although the models available to us are based on a carbon- 
rich chemistry, we nevertheless consider this comparison to be 
of some relevance, given that both types of objects, C-rich and 
O-rich dust-enshrouded stars with large IR excess, are located at 
the tip of the AGB. In particular, we are interested in whether 
the dimensions in the O-rich case compare with those of the C- 
rich case, although the details of the wind driving mechanisms 
are expected to be different. We should also note that the uncer- 
tain distances already give this comparison a very approximate 
character. 

The above wind models apply grey radiative transfer and 
consequently the resulting radii do not depend on wavelength. 
They vary in time, though, due to the stellar pulsation (pre- 
scribed by a piston approach for the inner boundary condition). 
Each wind model is characterized by a set of the stellar param- 
eters mass M, luminosity L, effective temperature T^, as well as 
by the carbon-to-oxygen ratio C/O of the atmospheric chemical 
composition and pulsation parameters period P and velocity am- 
plitude Av of the piston. Mass-loss rates have been determined 
by averaging typically over 20 periods. 

As representative cases we selected t wo models (w69 and 
w63, see row no. 52 and 23 of Table 3 in IWachteretalJl2002h 
showing mass-loss rates of 1.8 and 5.0 X lO'^Moyr"', respec- 
tively. The stellar parameters are masses of 1 .0 and 0.8 Mq, tem- 
peratures of 2800 and 2600 K, respectively, and a luminosity 
(both) of L=1OOOOL0. The other parameters are, in both cases, 
C/0=1.3, f=640d, and Av=5kms The observable CSE di- 
ameter of the first model (with the lower mass-loss rate) varies 
between approximately 2000 and 6000 solar radii over several 
pulsation cycles. The second model (with the higher mass-loss 
rate) has a diameter between about 5000/?© and 9000/?©. With 
values between 9QQQRq and 13OOO/?0 the observationally de- 
termined FWHM of the CSE of our OH/IR stars are larger than 
the CSE diameters of these models but generally of similar ex- 
tension. The model radii based on the adopted stellar luminos- 
ity and effective temperatures vary between 340-520/?© and be- 
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tween 400-6007?©, and are also smaller than our best- fit UD 
radii and smaller than those derived from the K - L calibra- 
tion by Lepine et al. (1995). These differences may be due to 
different stellar parameters and/or an overestimate of the stellar 
radius because of the presence of molecular layers as discussed 
in Sect. [3] 

/.From the 1612 MHz observations (te Lintel Hekke rt et alJ 

11991 ) we know the expansion velocities of the OH maser shells. 
Here, the value of Vexp = IT.Tkms-' for IRAS 13479-5436 
is larger than those of IRAS 14086-6907 and IRAS 17020- 
5254 (See Table [1]). Given that the stellar luminosities are 
thou ght to increase with increasing OH expansion velocities 
(e.g.r iSeve nster 2002) . IRAS 13479-5436 would have a higher 
luminosity than the other two sources, which is also consis- 
tent with its l argest K - L color index and the calibration by 
iLepine et all (Il995h . Also, iBaud & Habing (!l983) developed 
a relation between the main sequence mass and the OH out- 
flow velocity of log(M/M0) - (v'exp - 8)/16, suggesting the 
same main sequence mass of ~2 M© for IRAS 14086-6907 and 
IRAS 17020-5254 and a larger main sequence mass of ~4M0 
for IRAS 13479-5436 together with it s higher lum inosity. 

Additionally, the calibrations by iLepine et alJ d 19951) as a 
function of color index indicate an increasing optical depth 
of the CSE of OH/IR stars with increasing K - L color in- 
dex, which is supported by our observations (see Sect. O. A 
larger optical depth of the CSE indicates a lar ger integrated 
mass-l oss history of the star. The wind models bv lWachter et alj 
(I2002h . which are coupled to a stellar evolutionary code, in- 
dicate that the maximum mass-loss rate increases with initial 
mass, and that the shape becomes more spiky towards the end. 
Altogether we may speculate that the larger optical depth of the 
CSE of IRAS 13479-5436 compared to the other two sources 
may be a result of a higher initial mass and an evolutionary 
phase very close to the tip of the AGB evolution and thus 
the spike of the mass-loss rate. The different optical depths of 
the CSEs of IRAS 14086-6907 and IRAS 17020-5254, which 
have similar expansion velocities, may then rather be caused 
by an earlier evolutionary phase of IRAS 17020-5254 com- 
pared to IRAS 14086-6907 along the same initial mass track. 
Accordingly, IRAS 17020-5254 could even be an intermediate 
object, lying between a Mira-variable and a tip- AGB star. 

6. Conclusions 

Using the VLTI/AMBER instrument, we determined angu- 
lar diameters for the three OH/IR stars IRAS 13479-5436, 
IRAS 14086-6907 and IRAS 17020-5254 and their CSEs. Even 
though our visibility fits are based on a simple two-component 
geometrical model, they provide characteristic radial dimensions 
of the stars and their dust-rich shells, as well as constraints on the 
relative flux contributions from these different components. 

Our near-IR visibility data also indicate wavelength- 
dependent features, which resemble those detected in earUer 
AMBER observations of semi-regular and Mira-variable AGB 
stars, which are correlated with the positions of H2O and CO 
bands, and which are interpreted as the result of atmospheric 
molecular layers lying above the continuum-forming photo- 
spheres. This result confirms that atmospheric molecular layers 
are a common phenomenon among AGB stars of very different 
luminosities and mass-loss rates, alike. 

We also confirm that the circumstellar dust shells con- 
tribute significantly to the near-IR flux of our sources, unlike 
in the case of Mira variable AGB stars, which can well be de- 
scribed by dust-free model atmospheres at near-IR wavelengths 



(e.g. [Woodruff et al.ll2004l: iFedele et al.ll2005l: IWittkowski et all 
'2008*). This indicates a large optical depth of the CSE, which 
we interpret as the result of the "superwind" phase, the final 
10000 to 30000 years of AGB evolution, when the mass-loss 
rate increases by a factor of 10-100. This result also observation- 
ally confirms that oxygen-rich AGB stars develop a similar "su- 
perwind" as carbon-rich AGB stars, although their wind-driving 
mechanisms are not yet understood to the same detail. 

Even though the three stars have the same classification, in- 
terestingly they show somewhat different characteristics. This 
could be attributed to their different initial masses and slightly 
different stages of evolution. 
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